Correlated Supernova Systematics and Ground Based Surveys 
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Supernova distances provide a direct probe of cosmic acceleration, constraining dark energy. This 
leverage increases with survey redsliift depth at a rate bounded by the systematic uncertainties. We 
investigate the impact of a wavelengtli-dependent, global correlation model of systematics in com- 
parison to the standard local-redsliift correlation model. This can arise from subclass uncertainties 
as features in the supernova spectrum redshift out of the observer photometric filters or spectral 
range. We explore the impact of such a systematic on ground-based supernova surveys such as Dark 
Energy Survey and LSST, finding distinctive implications. Extending the wavelength sensitivity to 
1.05 /im through "extreme red" CCDs can improve the dark energy figure of merit by up to a factor 
2. 



I. INTRODUCTION 

Type la supernovae (SNe) are standardizable distance 
indicators that can be seen to high redshift, over at least 
70% of the lookback time of the universe. In a given 
galaxy, they occur on 10^-10'^ year time scales, so for 
10^° galaxies in the visible universe this implies one SN 
per second. A three-year deep survey could theoretically 
detect 10^ SNe so if the distance to each were known 
to 20% (0.46 mag; i.e. no standardization applied, just 
identification as Type la), then one could attain 0.01% 
distance measurements in 25 equal-error distance bins 
statistically. In this sense SNe are the most powerful 
cosmological probe statistically. 

However, as for any probe, the key issue comes down to 
systematic uncertainties not statistics. For SNe, such un- 
certainties include flux calibration, dust corrections, and 
population drift. Numerous articles have addressed one 
or another of these, for various experimental designs [ij . 
Some uncertainties can be removed by using spectropho- 
tometry rather than filter measurements; others can be 
ameliorated by observations spanning a long wavelength 
range; others by "like to like" comparison of time-series 
or spectral features 0, Q . 

To allow general comparison of different SN surveys, 
a phenomenological systematic model offers a compact 
way of summarizing the specifics of the measurement and 
analysis chain that delivers distances. It also provides a 
basis with which to identify previously unknown sources 
of uncertainty. Specifically, the model can contain pa- 
rameters describing the distribution of supernova sub- 
sets, e.g. magnitude offsets, that when fit can be checked 
for deviations or inconsistencies. 

A key feature of such a model would be a systemat- 
ics floor, reflecting that the large number of SNe men- 
tioned above should not in reality give an asymptotically 
small precision. This is fundamentally different from "self 
calibration" - adopting a fixed parametric form and si- 
multaneously fitting the error parameters along with the 
cosmology [l3|- Self calibration only holds if there is no 
uncertainty in the form, only the parameters. This is a 



weaker meaning of systematic uncertainty, and in gen- 
eral we do not have strong confidence in the form for 
population drift, for example. Therefore we consider a 
systematics model that explicitly includes a floor. 

A well- used systematics model was introduced by 
consisting of local correlations between SNe within red- 
shift bins of width 0.1. Large numbers of SNe in a bin 
reduce the distance error to a finite floor, but continued 
cosmological leverage could be obtained b y ad ding inde- 
pendent SNe, those in another redshift bin[l8|. The floor 
is given in magnitudes by 
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where d is the luminosity distance, z is the redshift of 
the bin center, and Zmax is the maximum redshift of the 
survey; we refer to this as the LH systematics model. An 
important feature is that the floor does not vanish even 
for z = 0. The cosmological parameter estimation results 
arc fairly insensitive to the bin width (scaling the error 
amplitude as the inverse square root of the width) within 
the range = 0.05 - 0.4 0. 

However, we want to explore beyond this local corre- 
lation model, where only SNe within a restricted redshift 
range are correlated with each other. One can imagine 
that some systematics might correlate all SNe, or all SNe 
above a certain redshift (e.g. if they rely on near-infrared 
observations). This can give a more global correlation, 
and we consider a model for that here. 

In Sec. [H] we present a motivated model for correlated 
systematic uncertainties, and in Sec. IIIII propagate its ef- 
fects through to dark energy cosmology fitting, especially 
the dark energy equation of state figure of merit. We in- 
vestigate in particular the influence of various survey and 
supernovae characteristics on the figure of merit. Sec- 
tion IIVI extends this to error fitting and to error models 
describing population drift. The results are summarized 
in Sec. IVl 
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II. CORRELATED SYSTEMATIC 

An error model with every supernova potentially cor- 
related with every other supernova is general, but one 
should use some physical motivation to guide the form 
of the correlations flOj. Our systematic model is based 
on SN spectral characteristics and says that errors, once 
they enter the measurements, affect all SNe above the 
entry redshift. We consider several different error con- 
tributions, each with their own entry redshift. This cor- 
relates that portion of the systematic uncertainty com- 
pletely among all higher redshift SNe. The error matrix 
can be viewed as being "laid down" in steps, with the 
first error source contributing af to all elements of the 
error matrix, the second error source adding (t| to the 
matrix block corresponding to higher redshift SNe, and 
so on for even higher redshift SNe. Effectively, the er- 
ror matrix has a "wedding cake" form, climbing in steps 
from the upper left to lower right. 

Such an error structure naturally arises as supernova 
light rcdshifts out of the range of observational sensi- 
tivity. For example, the SN la subclassification scheme 
of Branch et al. (2006) Q identifies four groups based 
on the equivalent widths of absorption features at rest- 
frame 5750 and 6100 A (from Si II), and the pres- 
ence of high-velocity Ca II IRS absorption at 8000 A. 
All three lines are available for z < 0.12 candidates 
from ground-based spectroscopic observations within the 
0.35 ^m < A < 0.9 /zm atmospheric window. 

The spectral features contribute to the standardiza- 
tion of the SNe as distance indicators; this can be viewed 
as calibrating the SN absolute magnitude M. As one 
by one the spectral features become inaccessible beyond 
z = 0.12, 0.47, and 0.57, supernovae at higher rcdshifts 
incur additional error floors from imperfect marginaliza- 
tion over each lost calibration parameter. This produces 
an error model where progressively higher redshift bins 
share additive correlated uncertainties. Note that even 
with measurements of these spectral features, the ensem- 
ble of supernova distances is still expected to have a cor- 
related error floor due to systematic limitations of such 
a classification scheme in general. 

The principles of a wavelength dependent, globally cor- 
related systematic can be extended to the larger set of 
spectral features cited in the literature as indicators of SN 
la diversity Moreover, since relevant features gener- 
ally have large equivalent widths (due to the high veloc- 
ities in the exploding supernova), they may show up in 
broadband photometry as well. The bounded observer 
window also induces redshift-dependent correlations in 
broadband photometric analysis from propagated uncer- 
tainties of SN and dust models at common restframe 
wavelengths. We therefore investigate application of such 
a correlated systcmatics model to surveys involving either 
broadband photometry or spectroscopic foUowup. 

Mathematically, the SN absolute magnitude is taken 
to be determined as a function of measurements of N 
spectral indicators. At higher rcdshifts, some of these. 



say N — i redder indicators, are unavailable. This leads 
to increased magnitude uncertainty, and a covariance in 
magnitude (and hence distance) errors of 
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where the intrinsic dispersion enters as a diagonal contri- 
bution crfj^t (generally dependent on the procedure used 
to determine M for each supernova), and 
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The indices a, b label individual SNe and Ha is the num- 
ber of spectral features observed for supernova a. The 
more spectral information, the fewer terms contribute to 
the sum and the smaller the overall uncertainty. However 
note that we always keep a base variance af. 

An example covariance matrix for three error sources 
(and five SNe in increasing redshift order) is shown below: 
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Another way to achieve the same result is to view the 
table of error contributions for each SN as a A^sn x -^err 
matrix. The product of this with its transpose will au- 
tomatically give a positive definite, symmetric matrix. 
That is, let 
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Then the product V = SS^ gives Eq. Here we see 
that each row gives the systematic error summary of a 
SN, and each column represents a different systematic er- 
ror contribution. So all SNe suffer from uncertainty type 
1. with magnitude ai (representing that even with the 
identified spectral features there will be residual, uniden- 
tified systematics); only the last three SNe (which could 
represent higher rcdshifts) suffer from type 2 uncertainty, 
and only the last SN (e.g. highest redshift) has type 3 un- 
certainty. Thus, the fourth SN gets contributions from 
type 1 and 2 uncertainties, while the fifth SN has sys- 
tematics of types 1, 2, and 3. This is easily generalized 
to A'orr contributions to the systematic uncertainty and 
to A^sN supernovae. Such an error list, arising from a 
linear combination of independent uncertainties (similar 
to a Karhunen-Loeve transform), produces the wedding 
cake type of covariance matrix seen in Eq. ([4]), and we 
call this the KL systematics model. 

The amplitudes of the uncertainties can be expressed 
either in terms of the individual cto or in terms of the 
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covariance matrix entries 

^-±-1 (6) 

1 

For definiteness, and to compare the results to the com- 
monly used, diagonal LH systematic model, we explore 
the case where 

^l = <ylui{z.)). (7) 

For the ground-based survey we assume (Tlh(-2^) = 
0.03 (l-|-z)/1.9 and here (zi) is the average redshift of the 
SNe within the range where ai (but not yet (Ti+i) enters. 
In the limit that iVcrr equals the number of independent 
error bins used in the LH model, the diagonal elements 
of the KL and LH covariance matrices are the same. 

III. EFFECT ON COSMOLOGY FIT 

We now explore the implications of the KL correlated 
systematic model on the cosmological parameter estima- 
tion. 
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A. Cosmology and Survey Parameters 

We consider a survey with 2000 SNe uniformly dis- 
tributed between z = 0.1 — Zmax, supplemented with an 
additional 300 SNe at z = 0.03 — 0.1. The values of Zmax 
considered are 0.6, 0.7, 0.8, 0.9. The distance modu- 
lus covariance matrix consists of a diagonal contribution 
from an intrinsic dispersion a-mt =0.1 for each supernova 
and a correlated part given by the KL systematic model. 
We explore several KL scenarios with different N^^t and 
entry redshifts for each error, but in all cases the ampli- 
tude of the systematics is fixed by Eq. The statisti- 
cal measurement uncertainty is taken to be subdominant. 
To speed up our analysis, we combine supernovae in 0.1 
redshift bins and use the average distance modulus with 
variance crf^^/rii, where rii is the number of SNe in each 
bin. (Such a procedure is found to make negligible dif- 
ference to the final result.) We include a CMB Planck 
prior on the distance to last scattering. 

The cosmological fit parameters of interest are the 
matter density flm in units of the critical density, and 
the dark energy equation of state parameters wq and Wa 
where w(a) = wq + Wa (1 — a). Spatial curvature is taken 
to vanish. The fiducial model is a flat ACDM cosmology 
with fim ~ 0.28. The combination Ai of the SN absolute 
magnitude and Hubble constant is also fit. We exam- 
ine the impact of different systematic error cases on the 
dark energy figure of merit F0M=1/ y^det Cwgwaj where 
CwQWa is the 2x2 submatrix of the parameter covariance 
matrix (inverse Fisher matrix). 

In Sec. IIV Al we carry out a fit rather than a Fisher 
matrix analysis, finding the results agree well. We also 
consider there the case where the ai are fit simultaneously 
with the cosmology. 



FIG. 1: The dark energy figure of merit (FOM) is plotted 
vs survey redshift depth Zmax for different cases of the global 
systematics model. Adding more error layers, each appearing 
above its corresponding redshift, lowers tiie FOM for surveys 
extending beyond those redshifts. The LH local systematics 
model has a worse FOM because SNe in different redshift bins 
do not see correlated contributions among its 9 errors. 



B. Impact of Survey Characteristics 

A first point of interest is the behavior of the cosmo- 
logical parameter estimation as we include various un- 
certainties, i.e. as A^err changes. When A'orr = 1, then all 
elements of the error matrix are identical, correspond- 
ing to completely correlated uncertainty across all SNe 
and all redshifts; this only affects the uncertainty on the 
absolute magnitude parameter Ai , not any of the cosmo- 
logical parameters. When the off-diagonal elements are 
not all equal, as for > 1, then the values of each di 
are relevant. 

The number of features contributing systematic uncer- 
tainties is important. Figure [T] shows the effect of both 
the number of error levels and the maximum depth of the 
survey on the FOM. 

For the Ncn = 2 case, we take the second error to 
come in at z — 0.1; this can either reflect the high- 
est wavelength (Ca II) spectral feature or be viewed as 
the difference between a separate, local (z < 0.1) SN 
survey that may be carried out spectrophotometrically 
(e.g. Nearby Supernova Factory [10[) and a high redshift 
ground based SN su rvey such as Dark Energy Survey 
(DES [13) or LSST Note that the A^err = 9 and 

LH systematics cases both have new errors appearing 
every 0.1 in redshift, and in fact have identical diagonal 
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terms in the correlation matrix. However, the global cor- 
relation causes more of the effect of the uncertainties to 
propagate into the absolute magnitude parameter, ame- 
liorating their influence on the cosmological parameter 
determination, and so giving a higher FOM. 

We take for our canonical systematic model the case 
mentioned in Sec. where the redshifting of important 
Si II and Ca II spectral features out of the observed wave- 
length range generates systematic uncertainties. For sen- 
sitivity diminishing at 0.9 ^m these occur at z = 0.12, 
0.47, 0.57. This is referred to as the A'on- = 4 case since 
the error matrix will be built up out of 4 layers, i.e. blocks 
involving through The following quantitative re- 
sults use this model. 

Most interestingly, the KL systematic, if it holds as a 
better description of SN observations for ground based 
observations than the uncorrelatcd LH systematic, pre- 
dicts (for iVoir = 4) a factor of 3 greater FOM for com- 
prehensive surveys extending to Zmax = 0.9. Moreover, 
the constraints in the dark energy equation of state wq- 
Wa plane have slightly more complementarity with other 
probes such as weak lensing or baryon acoustic oscilla- 
tions than in the LH case. 

Considering the number TV of well-characterized SNe 
in the survey, we find that the FOM does not saturate 
with number as severely as in the LH systematic case. 
Instead, we have 



FOM - N^/'^ . 



(8) 



over the range considered; for very large TV it gradually 
levels off. In the purely statistical error case FOM^ iV 
but recall that in addition to these N SNe at z = 0.1- 
^max we include a constant 300 SNe at z = 0.03-0.1 and a 
Planck CMB prior on the distance to last scattering. We 
emphasize that the iV SNe represent those with spectral 
or sufficiently high signal to noise photometric informa- 
tion to discern the different spectral feature classes down 
to the level of the Esys(z) = 0.03(1 -I- z)/1.9 systematic. 
Figure [5] shows the FOM as a function of N. 

Regarding the amplitude of the systcm- 
atics, in the canonical case, the errors are 
{crj = (0.0167,0.0117,0.0130,0.0133). The more 
correlated we make the SN uncertainties, by raising 
tTi, the less the effect of the remaining errors on the 
cosmological parameters. Figure [3] illustrates this, where 
we adjust cti to be either 0.02 (which drives (72 nearly 
to 0, i.e. S2 Si due to the constraint of Eq. [7]) or 
(making (T2 — 0.020). The area of the wo-Wa contour 
clearly decreases as more of the systematic becomes 
global, and interestingly the degeneracy direction rotates 
slightly to become more complementary to other probes 
such as weak lensing and baryon acoustic oscillations. 

Studying the highest redshift systematic uncertainty 
cr4, we find an analogous effect on the FOM, as seen in 
Fig. H As (74 ^- 0, S4 S3 and the SNe above z = 0.47 
all have a coherent systematic. This improves the FOM 
by 33%. Thus, tight control of systematics out to the 
highest redshifts is beneficial. As (T4 is instead increased. 
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FIG. 2: The dark energy figure of merit (FOM) is plotted 
vs number A'' of well-cliaracterized SNe in tire high-redshift 
survey. Tlie canonical FOM for = 2000 changes by factors 
1.5 or 0.7 for A'' doubling or halving, respectively. 
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FIG. 3: 68% confidence level contour in the dark energy equa- 
tion of state plane is shown for varying levels of the global 
systematic ai. The sum E2 = + (t| stays fixed. As 
CTi approaches E2 = 0.0204, 0-2 — and the SNe out to 
z = 0.47 share a common systematic, propagating mostly 
into A4 rather than dark energy parameters. Thus the con- 
tour tightens. Conversely, as ai — 0, the SNe become less 
correlated and the constraints weaken. 
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FIG. 4: The dark energy figure of merit (FOM) is plotted for 
different values of tfie highest redshift systematic (T4. If the 
spectral feature at 5750 A responsible for 0-4 is not actually 
relevant to SN magnitudes, so (74 — >■ 0, then the correlation in 
magnitude uncertainty among the highest redshift SNe grows 
and the FOM improves. Larger systematics lower the FOM 
at all redshifts beyond z — 0.6 where this error enters. Here 
af'^ is the canonical value 0.0133 given by Eq. ([7|. 



FIG. 5: The dark energy figure of merit (FOM) is plotted vs 
survey redshift depth Zmax for different values of the wave- 
length sensitivity cutoff. Increasing the wavelength range can 
give substantial improvements to FOM, even for shallower 
surveys. 



the FOM worsens, but slowly, with a doubling of (T4 (30% 
inerease in S4) causing an 8% decline in FOM. Since 
tT4 enters at z = 0.57, the FOMs for all values of 0-4 
agree up to this redshift; generally the FOM then flattens 
with survey depth as the impact of the new error enters. 
At even higher redshift the FOM gradually increases as 
the leverage from a larger redshift baseline asserts itself. 
Note however that the FOM can actually decrease with 
^max if the error is large enough to overwhelm the redshift 
leverage. 

Since the SN properties whose lack of measurement is 
responsible for the systematic uncertainties enter at ob- 
served wavelengths that increase for higher redshift SNe, 
an extension of the wavelength region over which the sur- 
vey is sensitive would forestall the loss of information due 
to missed features. This is of course one of the main rea- 
sons for considering space-based SNe surveys. But within 
the context of ground-based surveys, if wc could slightly 
extend the wavelength region over which high signal to 
noise measurements could be made, this could pay divi- 
dends for cosmological parameter reach. 

The ongoing development of "extreme red" CCDs p^ . 
detectors with enhanced efficiency out to l/xm or even 
1.05 ^m, offers the possibility of recovering from some of 
the observing inefficiencies caused by bright atmospheric 
emission and absorption at longer wavelengths. We con- 



sider these two cases in addition to the standard 0.9 /im 
cutoff. Figure [5] shows that such technology advances 
could have dramatic effects on dark energy science. 

The general effect of a systematic error contribution is 
to suppress the growth of FOM with redshift at the red- 
shift corresponding to where the spectral feature leaves 
the observational sensitivity window. That is, the curves 
in the figure flatten at the redshifts where the errors en- 
ter, as mentioned above. For the Amax = 0.9 /im case, 
all the errors have entered by z = 0.57 and so the FOM 
curve, although starting low due to the accumulated sys- 
tematics, increases with Zmax- For the Amax = 1 Mm case, 
the errors appear at z = 0.25, 0.64, and 0.74, hence the 
FOM curve is roughly flat over the range of z,„ax shown. 
For Amax = 1.05 /im, the errors come in at z = 0.31, 
0.72, 0.83 so the FOM is able to climb with redshift out 
to Zmax = 0.7 and then levels off. 

Compared to the canonical Zmax ~ 0.9 case the im- 
provement in FOM by extending sensitive measurements 
to Amax = 1 or 1.05 //m is 7% or 32% respectively. More- 
over a survey with extreme red sensitivity to 1.05 //m is 
equally powerful at Zmax = 0.7 vs 0.9 (recall the number 
of SNe is kept constant at = 2000), giving there a 
factor of 2 improvement in FOM over a standard 0.9 /im 
cutoff survey. 
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IV. FURTHER APPLICATIONS 
A. Error Fitting 

Systematics represent the residuals after known errors 
are taken into account. One usually estimates, either em- 
pirically through training and validation subsets or using 
simulation and modeling guidance, the level of the sys- 
tematic. An example is the 0.03 mag amplitude in the 
LH model. In Sec. IIIII we assume that those levels are 
already well-characterized when incorporated in the dis- 
tance modulus covariance matrix. Alternatively, within 
the cosmology fit one can simultaneously estimate the 
level of variance or introduce new parameters to fit out 
the effect of the systematic. 

To explore this we carry out simultaneous fits to vari- 
ous systematics model parameters along with the cosmo- 
logical parameters. In place of the Fisher analysis, we 
use MINUIT [ij] for parameter fitting and the resultant 
Hessian matrix to calculate the FOM. Contributions to 
the fit likelihood come from not only the usual term 
but also a penalty term In det C involving the covariance 
matrix C . For example, this approach was carried out 
to fit the intrinsic SN dispersion CTint by [l^. We do not 
bin but analyze the full uniformly-distributed 2000 SNc 
plus the 300 local SNc. Since we are interested in fit un- 
certainties, when possible we do not realize dispersion on 
the simulated data to ensure that the best fit and input 
parameters are identical. We restrict these studies to the 
= 4, Zmax = 0.9 case. 

We first check that this likelihood approach gives the 
same results as the Fisher analysis for our canonical case. 
The FOMs agree to better than 5%. 

Next, we consider that the estimate of the absolute 
magnitude may have a different intrinsic dispersion de- 
pending on the number of spectral features used in its 
estimate. This is done by considering the cr?jj „ in 
Eq. ^ as (four) fit parameters rather than assuming 
they equal 0.1. Since we are fitting for variances, we 
must instantiate dispersion in the simulated data. We 
perform only one instantiation as a single fit is time con- 
suming with a 2300 x 2300 matrix inversion in each like- 
lihood evaluation. The uncertainty found for each Uint.na 
is ^ few X 10~^, the level to which we can detect real 
differences in the absolute magnitude models. The FOM 
is within 3% of the fixed a-mt =0.1 case. 

The correlated systematics ct^ reflect the residual un- 
certainties in the absolute magnitude offsets in redshift 
bins defined by spectral-feature visibility. These A'cn- 
offsets can be fit for simultaneously with the cosmology 
since the luminosity distance model is smooth, without 
sudden steps (see the next subsection for the case of a 
smooth magnitude drift). Fitting for the biases on Ali 
(without the V term in the covariance matrix), the real- 
izations of 5Mi are determined to ~ 10~^ and the FOM 
falls to 59 (i.e. by ~ 40%). However, when adopting 
the ai given by Eq. (O as priors on 5Mi (with mean 0), 
the realizations of A'h arc determined to ^ 7.5 x 10~^ and 



the FOM improves to 102 (by 12%). The smooth Hubble 
diagram of the cosmological parameters allows the self- 
calibration of magnitude offsets, leading to an improved 
FOM over the canonical case. 



B. Population Drift 

As used so far, the KL model assumes that once in- 
formation on a spectral feature is lost, the uncertainty 
enters and remains constant with redshift. In terms of a 
magnitude correction based on that feature, this makes 
sense: no information means no correction and hence the 
added uncertainty. But we can also consider the spectral 
feature as aiding classification of subtypes of SNe la with 
different absolute magnitudes. 

The various subtypes will all combine to give some av- 
erage absolute magnitude for the sample. With the loss 
of information, the proportion of subtypes is unlikely to 
suddenly change abruptly, and so the average should ini- 
tially remain a good approximation. At higher redshifts, 
however, the population fractions may drift further and 
so the previous average becomes increasingly uncertain 
as an estimate of the true sample average. In this pop- 
ulation drift interpretation, the systematic uncertainty 
increases as the redshift increases beyond that where the 
feature is lost. 

We explore a modified KL model where the systematic 
scales as 



(9) 



in the error table S. The exception is cti which we take 
to be a global uncertainty for SN standardization. The 
covariance matrix then contains elements 



{6Ma 6Mb) = a\ 



N —m\n. (na,?ib) + l 



+ 



E 

i=2 



Zh 



(10) 



where Za^ Zb are the redshifts of the two SNe being cor- 
related and Zi is the entry redshift of the systematic ai. 

The systematic ai is thus not a constant entry in a 
wedding cake layer but rather gives a ramp up: it starts 
at when it enters at Zi and attains an amplitude cr^^ at 
z = Zi, = 0.9. For comparison we choose aii, to equal the 
ai as given for the redshift independent, or flat, KL model 
in the earlier part of the paper to match the LH model 
amplitudes. Due to the scaling with z, the systematic 
amplitude will be less than the flat KL case for all z < z^,. 
Not surprisingly, therefore the FOM for the scaling case 
is larger (by 28% in the canonical iVorr = 4, Zmax = 0.9 
case). This implies that such population drift systematics 
have less effect (for the same amplitude at the maximum 
redshift) than the magnitude correction systematic. 
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V. CONCLUSIONS 

Systematics are the ultimate limiting factor in any cos- 
mological probe. We show that uncertainties due to the 
loss of information of supernova spectral features can cor- 
relate distance errors for all redshift higher than where 
the information is lost. One expects such loss in ground- 
based surveys where sensitivity declines strongly beyond 
9000 A. 

Such correlation in the error model can actually have a 
tempering effect, relative to a localized uncertainty such 
as in the LH model, in that a wider range of supernovae 
experience the effect, it propagates mostly into the ab- 
solute magnitude determination, and so the impact on 
the cosmological parameter estimate is reduced. For the 
same level of dispersion (autocovariance) , the KL corre- 
lated systematic model allows a larger dark energy figure 
of merit and slower approach to the systematic floor. 

Extending the wavelength sensitivity into the near in- 
frared, such as may be possible with new technology "ex- 
treme red" CCDs that also cover the optical region, has 
science potential that should be considered in mission de- 
signs. For a survey depth of Zmax = 0.9 the gain in figure 
of merit is 32% on going out to 1.05 /J,m. Moreover, the 
same FOM can be achieved even if the survey depth is 
^max = 0.7, since the added systematics beyond z = 0.72 
due to loss of Si II features counteract the increased red- 
shift baseline, for a fixed number of supernovae. 

Viewing spectral information loss as opening the analy- 
sis to uncorrected population drift, we find that this does 
not have as much impact as for the magnitude correction 



ansatz. As redshift increases, information might also be 
gained from rest frame UV features, but this would occur 
at low redshift (z w 0.1-0.25) and is not expected to have 
much impact on survey depth design. Currently the sit- 
uation concerning UV features for magnitude correction, 
or their ability to substitute for the higher wavelength 
information discussed here, is not clear. 

The calculations here are meant illustratively, to show 
the method and possible impacts of correlated system- 
atics. We do not claim that the Branch et al. model Q 
for supernova subclassification is truth, that those par- 
ticular wavelengths carry the essential information on su- 
pernova magnitudes. For example. Bailey et al. have 
identified certain distinct series of spectral flux ratios and 
ranked their correlation with absolute magnitude. Never- 
theless, the approach of analyzing correlated systematics, 
and the benefit of extending the wavelength range, e.g. 
through red sensitive CCDs, should be general. 
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